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Abstract. We have carried out high-precision timing measurements of the binary millisecond pulsar PSR 
J2051--0827 with the Effelsberg 100-m radio telescope of the Max-Planck-Institut fiir Radioastronomie and 
with the Lovell 76-m radio telescope at Jodrell Bank. The 6.5-yrs radio timing measurements have revealed 
a significant secular variation of the projected semi-major axis of the pulsar at a rate of i = d{ai sin i) / dt = 
(—0.23 ± 0.03) X 10~^^, which is probably caused by the Newtonian spin-orbit coupling in this binary system 
leading to a precession of the orbital plane. The required misalignment of the spin and orbital angular momenta 
of the companion are evidence for an asymmetric supernova explosion. We have also confirmed that the orbital 
period is currently decreasing at a rate of Pb ~ (—15.5 ± 0.8) x 10~^^ s s~^ and have measured second and 
third orbital period derivatives d^Pb/dt^ = (-f2.1 ± 0.3) x 10"^° s~^ and d^P^/dt^ = (3.6 ± 0.6) x 10"^** s"^, 
which indicate a quasi-cyclic orbital period variation similar to those found in another eclipsing pulsar system, 
PSR B1957-(-20. The observed variation of the orbital parameters constrains the maximal value of the companion 
radius to _Rc max ~ 0.06 Rq and implies that the companion is underfilling its Roche lobe by 50 %. The derived 
variation in the quadrupole moment of the companion is probably caused by tidal dissipation similar to the 
mechanism proposed for PSR B1957-I-20. We conclude that the companion is at least partially non-degenerate, 
convective and magnetically active. 
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1. Introduction 

The eclipsing binary millisecond pulsar PSR J2051— 0827 
was discovered with the Parkes 64-m radio telescope 
in Australia in a 0.4 GHz survey of the southern sky 
(Stappers et al. 1996a). This system has one of the short- 
est known orbital periods, Pb — 2.4 hrs, and is moving 
in an almost circular, compact orbit, with a separation 
between the pulsar and its companion of only 1.03 Rq. 
Observations of the pulsar with the 64-m Parkes and the 
76-m Lovell radio telescope at Jodrell Bank at frequen- 
cies between 0.4 GHz and 2 GHz showed that at lower 
frequencies near 0.6 GHz the pulsar is eclipsed by the at- 
mosphere of the companion during approximately 10% of 
the orbital period, whereas at higher frequencies near 1.4 
GHz there are almost no visible eclipses in this system 
(Stappers et al. 1996a). High precision timing observations 
of PSR J2051-0827 revealed that the orbital period of the 
system is decreasing at a rate of Pb ~ —10^^^ s s^^ indi- 
cating a decay time of the system of only 25 Myr (Stappers 
et al. 1998). Early optical observations of the field of PSR 
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J2051— 0827 revealed that the amplitude of the compan- 
ion's light curve is about 1.2 mag and the companion 
is probably rotating synchronously around the pulsar so 
that one side is being heated by the impinging pulsar flux 
(Stappers et al. 1996b). The fit of photometry data to a 
model of a gravitationally distorted, low-mass secondary 
star that is irradiated by the impinging pulsar wind, 
has shown that the inclination of the system is greater 
than 30° and the maximum companion mass is 0.055 Mq 
(Stappers et al. 1999). Recent observations of the pulsar's 
companion using the Hubble Space Telescope (HST) Wide 
Field Planetary Camera have allowed to detect its "dark" 
side by Stappers et al. (2001b). Surprisingly, they detected 
a slight asymmetry in the companion's light curve indi- 
cating that a simple synchronously rotating companion 
is no longer a complete model. Fitting the same model as 
Stappers et al. (1999) indicated that more than 30 % of the 
pulsar spin-down energy is converted into optical emission 
and that the system is moderately inclined, i.e. i ~ 40°. 

Long-term timing observations of eclipsing (and hence 
often interacting) binary pulsars can reveal additional 
variations in orbital elements, which are important for 
understanding the evolution models of these systems and 
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of pulsars in general. The orbital period derivative may 
change with time in a quasi-periodic way similar to PSR 
B1957+20 (Arzoumanian et al. 1994) or may remain un- 
changed over longer time intervals like for PSR B1744-24A 
(Nice et al. 2000). 

We have carried out high precision timing observa- 
tions of the binary system PSR J2051— 0857 during the 
last 6.5 years. In addition to improving previously deter- 
mined timing parameters, we detect, for the first time, 
time derivatives of the orbital period of second and third 
order and the time derivative of the semi-major axis. We 
use these discovered orbital variations to constrain the size 
of the pulsar companion. We discuss possible additional 
relativistic and non-relativistic effects that could be the 
cause for the apparent variations of the orbit. 

2. Observations 

PSR J2051-0827 was observed with both the 100- 
m radiotelescope of the Max-Planck-Institut fiir 
Radioastronomie in Effelsberg and the 76-m Lovell 
radiotelescope at Jodrell Bank. Observations have been 
carried out since 1996 using the Effelsberg 100-m radio 
telescope soon after the discovery of the pulsar with the 
64-m Parkes radio telescope in Australia. Timing data 
have been acquired approximately once per month, with 
a few larger gaps due to unavailability of telescope time. 
We have collected about 350 individual times-of-arrival 
(TOAs) at center frequencies of 0.86, 1.4, 1.7 and 2.7 GHz. 
Most of the observations were performed at frequencies 
near 1.4 GHz. 

During a period in August-October 1996, the data were 
obtained using the Effelsberg Pulsar Observation System 
(EPOS). Here two channels of circular polarization are 
processed in a polarimeter and a 4 x 60 x 666 KHz filter 
bank combined with an incoherent hardware dedisperser 
(Kramer et al. 1997, 1998). Only the two orthogonal total 
power signals of 40 MHz bandwidth were used for later 
analysis. The input signal from the radio telescope is syn- 
chronously accumulated into 1024 pulse bins with the ap- 
parent pulsar spin period. The time stamp of each inte- 
gration is synchronised to time signals from a hydrogen 
maser clock calibrated with the signals from the Global 
Positioning System (GPS). The dedispersed pulsar pro- 
files were integrated for 15 s and later transformed to 8 
min - integrations for further template matching in the 
time domain, using the cross-correlation of the integrated 
profile with a high signal-to-noise template. 

Since October 1996, the data were collected using the 
Effelsberg - Berkeley Pulsar Processor (EBPP), a coher- 
ent dispersion removal processor, installed at Effelsberg 
(Backer et al. 1997). The EBPP provides 32 channels for 
both polarizations with a total bandwidth of up to 112 
MHz depending on the observing frequency and dispersion 
measure. For PSR J2051-0827 bandwidths of 28 MHz, 56 
MHz and 112 MHz were available at 0.9 GHz, 1.4 GHz and 
2.7 GHz respectively. The outputs of each channel are fed 
into dedisperser boards for coherent on-line dedispersion. 
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Fig. 1. Integrated profile of PSR J2051-0827, obtained 
with the EBPP at 1.4 GHz (left) and the profile template, 
constructed from fitting of composition of three gaussians 
to the integrated profile (right), which is used for template 
matching. 



All 64 output signals are integrated synchronously with 
the apparent pulsar period. The integrated profiles (5 min 
integration time) were matched in the frequency domain 
(Taylor 1991) with a synthetic template. The template it- 
self consists of three Gaussian components (Fig. 1 right) 
that were fitted to a high signal-to-noise pulsar profile 
(Fig. 1 left) using the method developed by Kramer et al. 
(1994, 1998). The time stamp of each integrated profile 
was provided by a hydrogen maser clock calibrated with 
signals from the GPS to the Universal Coordinated Time 
UTC(NIST). 

In order to get a better timing solution for the orbital 
parameters of the system, the sets of TOAs in several ob- 
serving sessions were acquired during one complete pulsar 
revolution around its orbit, equal to ~ 2.4 hrs. Typical 
TOA uncertainties in observations are ~ 8 /iS at 1.4 and 
1.7 GHz, - 10 ^is at 2.7 GHz and - 25 ^is at 0.9 GHz. 

The Jodrell Bank data were collected since August 
1994 using cryogenic receivers at 0.3, 0.4, 0.6 and 1.4 GHz. 
The observations were made using both circular polariza- 
tion directions utilising a 2 x 64 x 125 KHz filter bank. 
The added output signals from both polarizations were 
filtered and digitised for further on-line hardware dedis- 
persion and synchronous folding with the pulsar's appar- 
ent period. Integrated profiles with integration times of 
3 minutes were matched by standard pulse template to 
obtain topocentric TOAs (BeU et.al 1997). Typical TOA 
uncertainties in these observations are ~ 10 /is at 0.4 and 
1.4 GHz and - 16 /iS at 0.3 and 0.6 GHz. 



3. Data reduction and Timing solution 

The obtained TOAs corrected to UTC(NIST) were fit- 
ted to a spin-down model of the pulsar rotation in the 
binary system using both the TIMAPR [] (Doroshenko 
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Table 1. Spin, astrometric, Keplerian and post-Keplerian 
parameters of binary pulsar system PSR J2051— 0827. 




1994 1995 1996 1997 1998 1999 2000 2001 
Year 



Fig. 2. Post-fit TOA residuals of PSR J2051-0827 plot- 
ted as a function of epoch, after subtracting a global 16- 
parameter timing model from the data. 

& Kopeikin 1995) and the TEMPO software package, 
which both use the DE200 ephemeris of the Jet Propulsion 
Laboratory (Standish 1990) and various relativistic tim- 
ing models. Since the pulsar is moving in an almost cir- 
cular orbit, we have applied the Blandford and Teukolsky 
phenomenological timing model (Blandford and Teukolsky 
1976). We have also fitted the model of Laplace-Lagrange 
parameters for systems with small eccentricity (Lange et 
al. 2001) using TEMPO. Comparing results obtained by 
TIMAPR and TEMPO are in excellent agreement within 
their uncertainties at a Icr-level. 

In the fitting procedure the three TOA segments, ob- 
tained with both the Effelsberg and Jodrell Bank data 
acquisition systems, were fitted for a mutual offset. Data 
lying in the range of orbital phases between 0.20 and 0.35 
were excluded from the fit because of the additional mod- 
ulation of the excess column density on the TOAs in the 
eclipsing region (Stappers et al. 1998, 2001a). From the fit 
we obtained precise estimates of the pulsar's parameters 
and those of the orbit. These include astrometrical param- 
eters, spin parameters, and Keplerian and post-Keplerian 
orbital parameters (Table 1). 

Post-fit residuals are plotted as a function of observ- 
ing epoch (Fig. 2) and as a function of orbital phase (Fig. 
3). The upper limits on the values of DM, i, e and P, 
presented in Table 1, were obtained by the individual in- 
clusion of the corresponding parameter in the fitting pro- 
cedure. These parameters, as well as the longitude of pe- 
riastron wo, were fixed at zero while fitting for the other 
parameters. 

Since the proper motion in right ascension and dec- 
lination {fJ-aj IJ-s) has significant mutual covariances and 
the TOA uncertainties are relatively large, the global fit 
failed to give good estimates and uncertainties for fia and 
jjs- To obtain more precise values of the proper motion 



Parameter 



Value 



Astrometric, spin and medium 
Right ascension a (J2000) 
Declination S (J2000) 
Proper motion in R.A. fia (masyr^^) 
Proper motion in Decl. fis (masyr~^) 
Spin period P (ms) 
Period derivative P (s s~^) 
Spin frequency u (s~^) 
Spin frequency derivative u {s~^) — 
Dispersion measure DM (pccm"'^) 
Derivative of DM, DM (pc cm"^ yr~^) 
Epoch (MJD) 



parameters 
20''51"07''.514 5(1) 
08°27'37".795(5) 
5.3(1.0) " 
0.3(3.0) " 

4.508 641 744 971 6(2) 
1.2737(5) X 10"^° 
221.796 287 344 250(8) 
6.266(2) X lO-^** 
20.7449(4) 
0.001 1(3) 
51 000.0 



Parameters of binary 
Orbital period Ph (d) 
Projected semi-major axis x (It-s) 
Eccentricity e 

Epoch of periastron To, (MJD) 
Longitude of periastron, uq 
Orbital period derivative Pb (s s~^) 
Second derivative of Pb, Ph (s~^) 
Third derivative of Pb, d^P^/dt^ (s"^) 
Derivative of semi-major axis i (s s^^) 



orbit 

0.099 110 250 6(2) 

0. 045 052(2) 

< 0.000 096 

50 999.983 601 7(9) 

0°.0 ' 
-15.5(8) X 10"" 

2.1(3) X 10"^" 

3.6(6) X 10"^** 
-23(3) X 10"" 



Upper limits to parameters 
Second derivative of x, x, (s"^) 7(15) x 10~^^ 

Derivative of eccentricity |e| (s"^) 1 x lO"^'^ 

Second Period deriv. P (s"^) 5(2) x 10"^^ 

Second frequency derivative v (s"'') 2(1) x 10"'^'' 

Deduced and additional parameters 



Distance (DM - based), (kpc) 
Transverse velocity Vt, (km s~ 
Shklovskii effect Pt (s s"^) 
Mass function fm (Mq) 
Companion mass nic (Mq) 

rms TOA residual a (ps) 
Timing data span (MJD) 
Number of TOAs 



1.3 

33(19) 

0.01 H- 0.09 x 10"^° 
0.999 50(1) X 10-5 
0.0273/ sin i " 

21.1 

49573-51908 
584 



" See Fig. 4 for the values and uncertainties in proper motion. 

These parameters were fixed at zero in global fit. 

Calculated assuming a pulsar mass rup = 1.4 Mq. 
Note — Values in parenthesis are la uncertainties in the last 
quoted digit. 



parameters, we have mapped out the region of space 
near its global minimum in the (/ia, ^s) plane. Contours of 



Xq , where Xo is the minimum 
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of the statistic x^j are plotted in Fig. 4. They show regions 
of Icr (Ax^ = 2.3), 2cr (Ax^ = 6.2) and 3cr (Ax^ = 11.8), 
which correspond to 68.3%, 95.4% and 99.7% confidence 
for /ic and ^g. Within Icr - uncertainty, the proper motion 
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Fig. 3. Post-fit TOA residuals of PSR J2051-0827 plot- 
ted as a function of orbital phase of the pulsar in the bi- 
nary system with Pb — 2.4 hrs, after subtracting a global 
16-parameter timing model from the data. The TOAs ly- 
ing in the eclipsing region with the orbital phase (0.20- 
0.35) were omitted from the global fit. 

in right ascension is /Zq, = 5.3 ± 1.0 mas yr~^ and in decli- 
nation is = 0.3 ± 3.0 mas yr~^. The composite proper 
motion of PSR J2051-0827 is = 5.3 ± 3.0 mas yr^^. 




2 3 4 5 6 7 8 
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Fig. 4. Contours of near its global minimum in the 
(liajlJ's) plane. The contours enclose regions of 68.3%, 
95.4%, and 99.7% confidence for values of fia and /i^. The 
dot shows a low confidence value, corresponding to the 
minimal value of ■ 

The newly obtained celestial coordinates and binary 
parameters for PSR J2051— 0827 are in excellent agree- 
ment to previously published values. The long time span 
even allows to detect higher order time derivatives of the 
orbital period Pb and the projected semi-major axis of the 
pulsar X = Op sin i. The measured values of Pb, -Pb and i 
are highly significant and were detected at a level of about 
8(7, reducing the post-fit TOA residuals by ~ 25%. 



4. Variation of the pulsar orbit 

4.1. Possible effects of DM variations 

Variations of the dispersion measure, DM, which can be 
large in eclipsing binary systems, can have a significant 
influence on the TOAs and thus on the fitted orbital pa- 
rameters and their time derivatives. Therefore we test if 
the observed variations of the orbital parameters could be 
induced by DM variations. 

-3 I ' . ' . ' . ' . ' 1 

-4 - 
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-6 
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Fig. 5. Secular variations of the projected semi-major 
axis X = Op sin i of PSR J2051— 0827, obtained from fits 
of the data with excluded corresponding orbital phases. 
The solid line represents x — —7.3 /isyr^^ from Table 1 
and the dotted lines show the uncertainty in x. 

The best fit solution for PSR J2051-0827 reveals 
a secular variation of the dispersion measure at a rate 
d(DM)/dt = 0.0011(1) pc cm-3 yf-i (Table 1). This value 
is consistent with the expected square root dependence of 
DM variations to the distance, caused by the motion of 
the pulsar through the interstellar medium (Backer et al. 
1993). Apart from this long-term variation in DM, there 
might be variations over the orbit, which are due to a 
non-isotropic distribution of the electron column density 
rio near the pulsar eclipsing region. The observed value of i: 
may be caused by changes of ric near regions of upper and 
lower conjunction. Such variations depend on the orbital 
phase. If the assumption about the variations in rio near 
upper and lower conjunction is correct, the fit to the tim- 
ing model using the TOAs with excluded data segments in 
the regions of upper and lower conjunction, should show 
no variation in the projected semi-major axis x. Fig. 5 
shows these sets of timing solutions for i, where we have 
excluded 20% of the data around the corresponding or- 
bital phases (e.g., in obtaining solution of x at orbital 
phase 0.3 we have excluded all TOAs, lying in range of 
phases 0.2-0.4). It is seen that the detected secular change 
in X does not significantly depend onto the distribution of 
the electron column density along the orbit, so that DM 
variations cannot explain the observed x. 
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Date 

Fig. 6. (a) Fractional changes of the orbital period in 

the PSR J2051-0827 system vs. date, obtained from fit 
of the pulsar parameters within four segments of data. 
(6) Fractional changes of the projected semi- major axis 
of PSR J2051-0827 vs. date, (c) Fractional pulse period 
changes vs. date. The solid lines corresponds to the values 
of Pb, Pb and X listed in Table 1. 



We have also performed a more general test of the na- 
ture of orbital parameter variations in a way described by 
Arzoumanian et al (1994). The total data set was divided 
into five sub-intervals, each spanning about one year of 
TOAs. The individual TOA sub-sets were fitted for the 
pulsar spin and Keplerian parameters with variations of 
orbital parameters held fixed at zero. The time of peri- 
astron To and the orbital period Pb were transformed to 
and held fixed at an epoch near the center of each sub- 
set. The values of the pulsar spin period were transformed 
to the same epoch using the spin-down model of pulsar 
rotation, so that variations in period were calculated as 
APo = AP — PAt. The resulting fractional changes in 
the orbital period, the projected semi-major axis, and 
the spin period are shown in Figs. 6a-6c. Variations in 
DM should affect measurements of the spin period Pq, or- 
bital period Pb and semi-major axis x in the same way as 
APb/Pb = Ax/x = APq/Pq. As relative changes differ by 
a few orders of magnitude, we conclude that the variations 
in DM do not affect variations in Pb, x and P. 



4.2. Origin of tine variation of the orbital period 

We have confirmed the existence of a large orbital period 
derivative of the PSR J2051— 0827 binary system, previ- 
ously reported by Stappers et al. (1998). The observed 
vahie is equal to Pb = (-15.8 ± 0.3) x lO^^^ gy-i (g^^; 
Table 1) and is in good agreement with that found by 
Stappers et al. (1998). 

There are a number of effects that may cause the 
change of the orbital period of a binary system, which 
can be summarised as 



(1) 



where the contributions to the observed orbital period 
derivative are due to the emission of gravitational radi- 
ation {GW), acceleration of the binary system (ace), mass 
loss from the system (m), tidal dissipation of the orbit 
(r), and gravitational quadrupole coupling (Q). 

The change in orbital period due to the general rela- 
tivistic orbital decay, 




pGW =_(i95V5)(2^r0/Pb)5/3(l + 

-7/2 



73p2 
24 



37 4\ 
96^ J 



(2) 



X (l — e^) {mpmc)/{mp + nic)^/^ 



where Tq = GMq/c^ = 4.925490947 x 10"^ s (Damour & 
Deruelle 1986, Damour & Taylor 1992) and a pulsar mass 
of TTip = 1.4 Mq , is expected to be in the range = 
(—6.9 — 3.4) X 10~^^ s s~^ for the allowed companion 
masses rric = (0.027-^ 0.055) Mq . This value is about two 
orders of magnitude less than the observed value of Pt,. 

An acceleration of the binary system with respect to 
the Solar System Barycentre (SSB) may be caused by the 
differential rotation of the Galaxy or by a third massive 
body in the vicinity of the binary system. The acceleration 
affects the binary period derivative Pb and the spin period 
derivative P in the same way. Even under the assumption 
that the observed value of the spin period derivative is 
totally due to acceleration, the contribution to the binary 
period derivative will be only (Pb/Pb)"'='= = {P / PY"" = 
3 X 10~^* s~^, which is four orders of magnitude less than 
the observed value. 

The mass loss from the binary system leads to a change 
in the orbital period at a rate (Damour & Taylor 1991) 




Stt^ 



L. P.. 



mp + TUc (? Pi ' 



(3) 



where denotes the moment of inertia of body a (a = p 

for the pulsar, a = c for the companion). Pa and Pa are 
the spin period and period derivative. Since the compan- 
ion is almost synchronously rotating around the pulsar 
(Pc ~ Pb), we can neglect its contribution. Using observed 
values of the pulsar spin period and its time derivative, 
and assuming Ip ~ 10*^ g cm~^, the change of the orbital 



period due to mass loss is P™ 
hence too small to be of importance. 



4 xlO"^'^ s s~^ and 
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Tidal torques cannot be the reason for the observed Pb 
either, because the magnitude of the tidal torque is orders 
of magnitudes too small to transfer the necessary angular 
momentum (Applegate 1992, and references therein). 

As all other contributions are several orders of magni- 
tude smaller than the observed (Pb/Pb), we conclude that 
the change of the orbital period is most likely due to grav- 
itational quadrupole coupling which has been proposed 
earlier for the eclipsing binary system PSR B1957-f20 
(Applegate and Shaham 1994). A variable quadrupole mo- 
ment which is due to a cyclic spin-up and spin-down of 
the outer layers of the companion, provides a natural 
explanation of the quasi-cyclic orbital variations of the 
PSR B1957-I-20 binary system found by Arzoumanian et 
al. (1994) and recently confirmed by Nice et al. (2000). 
In fact, our observations revealed that the orbital pe- 
riod derivative of PSR J2051— 0827 changes with time 
at a rate of Pb = (+2.1 ± 0.3) x lO^^o g-i rpj^jg 
plies that the presently decreasing orbital period will in- 
crease after At = — 2Pb/Pb ^ 33 yrs and that the system 
is undergoing quasi-cyclic variations like that found for 
PSR B1957-I-20. The orbital period change corresponding 
to a variation of the quadrupole moment of the compan- 
ion, AQ, is equal to (Applegate and Shaham 1994) 



APb\^ 
) 



-9 



AQ 



(4) 



where rric is the companion mass, and a is the relative 
semi-major axis of the orbit (a = ap + ac). The semi- 
major axes of the pulsar ap = x/sini and companion 
Oc = a;(l. 4/0. 0273)/ s'mi are calculated for the pulsar and 
companion masses rup = 1.4 Mq and tUc = 0.0273 AIq 
respectively. The quadrupole moment of the companion 
can be derived as (Kopal 1978) 



9^ 



(5) 



where k is the apsidal motion constant (Claret & Gimenez 
1991), fie is the angular velocity of the companion, its 
radius, and G the Newtonian gravitation constant. By dif- 
ferentiating Equation (|^) the variation of the quadrupole 
moment with Qc is 



AQ 



4 0|i?| /Aa 
9 G \ Qc 



(6) 



From Equations (Q) and (||) the variation of the compan- 
ion angular velocity corresponding to the variation in the 
orbital period is derived as 



'n7 



1 / APb 
4 ^ Pb 



G 



rrica 



(7) 



We estimate that the radius of the pulsar's companion 
is Rc ~ 0.06 i?0 (see Section O). This value is equal to 
half the size of its Roche lobe Pl — 0.13 P0 (Stappers et 
al. 1996a) and agrees with that inferred from optical ob- 
servations of the companion (Stappers et al. 1999, 2000). 



The apsidal motion constant, fc, strongly depends on 
the effective temperature of the companion. The first 
optical observations indicated a temperature of Pcff — 
4000 ~ 4700 K (Stappers et al. 1996b) whereas recent ob- 
servations with the HST indicated that the backside tem- 
perature is likely to be less then 3000 K (Stappers et al. 
2000, 2001b). The corresponding value of the apsidal mo- 
tion constant is fc = 0.044^0.159 (Claret 2000). Since the 
companion spin frequency is close to that of the orbit, we 
can write for the general case 



27r 



(8) 



where / is close to unity. While / = 1 describes syn- 
chronous rotation, the wind and the magnetic activity of 
the companion in the Applegate & Shaham (1994) model 
result in a torque which slows down companion's rotation, 
so that / < 1. From the light curve observed by Stappers 
et al. (2001b, see their Fig. 1), we can expect this devi- 
ation from synchronous rotation to be small, estimating 
a lower limit of / > 0.9. For / « 1, the change of the 
orbital period of APb ^ —1.5 x 10^* s observed after six 
years is caused by a variation of the angular velocity of 
the companion of 



A17c 



(8 + 60) X 10' 



(9) 



for an orbital inclination 30° < i < 90° and a range of 
companion masses of nic ~ 0.027 + 0.055^/©. For / « 0.9 
we derive Afic/fic (1 + 74) x 10~^. These estima- 
tions of the angular velocity variations agree well with 
ArJc/f^c ^ 10"^ for the companion of PSR B1957+20 
and with Af7c/rJc - (3 + 30) x lO'^ for Algol, RS CVn 
and CV systems (Applegate & Shaham 1994). 

4.3. Origin of the variation of the projected semi-major 
axis 

We have detected a significant secular variation of the 
pulsar projected semi-major axis at a rate of i = 
d{asini)/dt = (-23 ± 3) x lO^^"* ss^^ in a global fit of 
the pulsar parameters. For an independent check of this 
effect, the TOA data set was split in five approximately 
equal segments, which were fitted for pulsar astrometric, 
spin and orbital parameters with the time derivatives of 
the orbital elements held fixed at zero. In these segments 
the epochs of ascending node were converted to values cor- 
responding to the centre of each data set. The fitted values 
of the projected semi-major axis for these three time in- 
tervals are found to be in good agreement with the timing 
solution obtained from the global fit, and confirm the cur- 
rent decrease of the projected semi- major axis (Fig. 7). 

The value of x cannot be explained by the observed 
time derivative in the orbital period Pb. The observed 
value of Pb assumes a corresponding change in the pro- 
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Fig. 7. Secular variation of the projected semi-major axis 
X = apsini of PSR J2051-0827. The dots indicate the 
values with their formal errors, obtained from fits of sub- 
sequent TOA data segments with the time derivatives of 
the orbital elements held fixed at zero. The solid line shows 
the secular variation of x obtained from the global fit. 



jected semi-major axis as a consequence of Kepler's third 
law, at a rate 



X 
X 



2I\ 
3Pb' 



(10) 



For the measured values of x, Ph and Pb one obtains 
imax — —6 X 10~^^s s~^, which is about 4 orders of mag- 
nitude less than the observed x. 

There are a few effects which may cause a secular 
change of the projected semi-major axis x (Damour & 
Taylor 1992, Kopeikin 1994, Wex & Kopeikin 1999), i.e. 



obs 



GW 



PM , 
X \ dEA 



dt 



D 
D 



so 



,(11) 



where the contributions to the observed semi-major 
axis derivative are due to the emission of gravitational 
radiation (GW), proper motion of the binary system 
(PM), varying aberration {dsA/dt), changing Doppler 
shift {~D/D) and spin-orbit coupling [SO) in the binary 
system. 

The first contribution (GW) in Equation ( pi] ) repre- 
sents the shrinking of the pulsar orbit due to gravitational- 
wave damping in the binary system given by 



GW 




GW 



(12) 



With P^^ calculated from Equation (||), one expects a 
rate of {d^/a^)'^^ - (-5.4 ^ -2.6) x 10'^^ s^^, which 
is about 6 orders of magnitude smaller than the observed 
(x/x) and can be neglected. 



The second contribution {PM) in Equation (11) is 
caused by the proper motion of the binary system 



(Kopeikin 1996, Arzoumanian et al. 1996; Bell et al. 1997) 
given by 



PM 



1.54 X 10"^^ coti ( — /Xq, sin -I- /ia cos ri) , (13) 



where $7 is the longitude of the ascending node of the or- 
bit, and ^.Q. and fis (mas yr^^) are the proper motions in 
right ascension and declination. In this case, the observed 
change in x is caused by the secular variation in the incli- 
nation angle i due to the proper motion, whereas the in- 
trinsic semi-major axis of the pulsar Cp remains constant. 
The maximal contribution of the proper motion is 



PM 



1.54 X 10"^^ M coti 



2^1/2 



(14) 



5(3) 



with a composite proper motion = (^^ -|- fif) 
mas yr~^. In order to estimate this quantity we still need 
an estimation of the maximal inclination angle i. 

If the companion is a white dwarf, we obtain a min- 
imal radius using the mass-radius relation of Paczyhski 
(1967) of i?cmin — 0.03 Rq for mass ranges of TOc — 
0.027 0.055 Mq. If the companion is, instead, a non- 
degenerate star in the pulsar radiation field, as suggested 
by Ergma et al. (1998), the radius is even larger. We hence 
obtain i?c min — 0.03 Rq as a safe lower limit. Since at 
higher frequencies the pulsar signal is not eclipsed by the 
companion that is separated by 1.03 Rq, we derive an 
upper limit for the inclination angle of imax — 88°. This 
inclination angle leads to i^ax — 1-2(7) x lO"^'* s s^^. For 
the likely inclination oi i ^ 40° (Stappers et al. 1999) we 
obtain x^^ = 4(2) x lO^^"^ s s^^. This is at least four 
orders of magnitude smaller than the observed x and can 
be neglected. 

The third term, (dsA/dt), in Equation ( pT|) is due to 
varying aberration caused by relativistic precession of the 
pulsar spin axis 



dSA 

dt 



P 

K (T 



Qgeod 



[cotX sin 2?7 + cot i cos rj) 



2)1/2 



sin A 



(15) 



where P is the pulsar spin period, and (?;, A) are the po- 
lar coordinates of the pulsar spin (see Damour & Taylor 
1992). The rate of geodetic precession given by (Barker 
and O'Connel 1975) 



(-12/3 
Qgeod _ ^ 

P "2c2(l-f 



5/3 



mc(4mp -I- 3mc) 
(TOp -t- Tnc)4/3 



(16) 



amounts to fis""'^ - 1.7 x IQ-i" = 0.3° yr^^ for a com- 
panion mass of nic = 0.055 Mq. This leads to (dsA/dt) < 
8.5 X 10~^° s~^, which is about 6 orders of magnitude less 
than the observed (x/x) and can be neglected. 

The fourth term in Equation (ITTI), (D/D), is due to a 
changing Doppler shift (D) caused by the change of the 
distance between the SSB and the binary pulsar system 
due to the acceleration of the binary system in the gravita- 
tional field of the Galaxy and due to the proper motion on 
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the plane of the sky (Shklovskii 1970, Damour & Deruelle 
1986). We obtain 



D 
D 



Ko • (apsR — agse) 



(17) 



where Kq is the unit vector of the pulsar along the line 
of sight, apsji and agge are the Galactic accelerations at 
the location of the PSR J2051— 0827 binary system and 
the SSB, Vt the transverse velocity of the pulsar, and 
d the distance between the pulsar and the SSB. Using an 
estimation of assB 2.1 x lO^^m s^^ (Carlberg & Innanen 
1987) we derive for the first term in Equation (|l^) an 
upper limit of 7 x 10"^^ s^^. The transverse velocity of 
the pulsar, Vt — 33kms^^, and a distance to the pulsar 
of d = 1.3 kpc yield a value of y|/cd - 10" s"^ Thus, 
the total contribution from a varying Doppler shift to the 
observed {x/x) is less than 0.3% and can be neglected. 

As all other contributions are several orders of magni- 
tude smaller than the observed [x/x), the variation of x 
is most likely caused by classical spin-orbit coupling (SO) 
in the binary system (Equation ^) . Spin-orbit coupling 
leads to a variation in the inclination angle, (di/ dt), while 
the semi-major axis Op remains constant. In this case the 
observed value [di/dt) is 



di 
dt 



so 



tani 



(18) 



If we assume that the observed x is totally due to spin- 
orbit coupling, we obtain (di/dt) = —0.008° yr"^ for an 
orbital inchnation of i = 40° and (di/dt) = -0.26° yr^^ 
for i = 88°. 

The rotationally induced quadrupole of the companion 
produces a variation of the orbital inclination angle as (Lai 
et al. 1995, Wex 1998) 



di 1 ^ kRlVil . . ^ 

dt = 2''-' aHl-e^)^ sm(2g)sm<I>. 



(19) 



Here fiorb = 27r/Pb is the orbital frequency and Qs = 
flc/(Gmc/ Rc)^^^ the dimensionless spin of the compan- 
ion, Rc its radius, a the semi-major axis of the relative or- 
bit, k the apsidal motion constant, 9 the angle between the 
companion spin angular momentum S and the orbital an- 
gular momentum L, and $ the orbital plane precessional 
phase. Equation ( [l9| ) is valid for the case that |S| <C |L|, 
which is true for PSR J2051-0827, where |S| ~ 10-3|L|. 
Note that spin-orbit coupling requires the companions's 
spin axis to be inclined with respect to the orbital angu- 
lar momentum vector (9 ^ 0). Measurement of (di/dt) 
and the classical periastron advance lu due to spin-orbit 
coupling can in principle be used to obtain constraints 
of the values of 9 and $ (Kaspi et al. 1996, Wex 1998). 
Unfortunately, the PSR J2051-0827 system has a negli- 
gible small eccentricity and all effects of secular variations 
in oj, along with the substantial relativistic advance of 
periastron oj^^ ~ 12° yr"-'^, are fully absorbed by the re- 
definition of the binary period. Therefore, they are not 
observable in this system (Kopeikin & Ozernoy 1999), so 



that we cannot obtain constraints from measurements of 
LU. Nevertheless we can yield an upper limit for (di/dt) 
and thus for the companion's radius from equations (iq). 



tan 1,1 



(20) 



where the indices max, min refer to the maximal or min- 
imal value of the corresponding parameter substituting 
observed values of imax and Xmin from Table 1. Using the 
maximal inclination imax = 88°, the corresponding upper 



limit (di/dt)^ 



-1.7 X 10-1" s-i, k = 0.044 0.159, 



a ~ 1.03 Rqi and e = 0, we then obtain for the case of 
a synchronously rotating companion (fJc — f^orb) a max- 
imal companion radius of i?cmax ^ 0.05 Rq. According 
to optical observations, the system is likely to be moder- 
ately inclined with an inclination angle i ^ 40° (Stappers 
et al. 2000, 2001b) which would lead to (di/dt)^^^ 
—4.8 X 10~i^ s~i and a maximal companion radius of 
Rc max ^ 0.03 Rq. The latter value is even smaller than 
the miminal radius for a pure helium white dwarf com- 
panion (Paczynski 1967) and therefore not very likely. 
If the companion is not perfectly synchronously rotating 
(/ ~ 0.9, see Equation^), we obtain -Rcmax ^ 0.04 i?Q for 
i - 40° and i?cmax ~ 0.06 Rq for i ~ 88°. A radius of 
Rc ~ 0.06 Rq is equal to half the size of its Roche lobe 
i?L = O.13i?0 (Stappers et al. 1996a) and agrees well with 
that inferred from optical observations of the companion 
(Stappers et al. 1999, 2000). 

5. Discussion and Conclusions 

The detected quasi-cyclic orbital variations of the 
PSR J2051— 0827 binary system are most likely caused 
by a variable quadrupole moment of the companion. For 
PSR B1957+20 Applegate & Shaham (1994) suggested 
that the companion's wind due to the pulsar irradiation 
and its magnetic activity result in a torque that forces 
the companion slightly out of synchronous rotation. We 
propose that this mechanism is the likely source of orbital 
variability also in the PSR J2051— 0827 binary system. In 
this framework the resulting tidal dissipation of energy is 
the source of the magnetic activity that causes the cyclic 
spin-up and spin-down of the outer layers of the compan- 
ion. This makes the companion of PSR J2051— 0827 the 
only second identified tidally powered star. We conclude 
that the companion is at least partially non-degenerate, 
convective and magnetically active (cf. Arzoumanian et al. 
1994) in accordance with a model proposed by D'Antona 
& Ergma (1993). 

The torque acting on the companion in the Applegate 
& Shaham (1994) model causes small deviations from a 
perfect co- rotation (/ < 1). This situation may be an 
explanation for the slight asymmetry in the companion's 
light curve that was recently observed by Stappers et al. 
(2001b). We note, however, that even though a similar 
asymmetry detected for the companion of PSR B1957-I-20 
by Djorgovski & Evans (1988) could not be confirmed by 
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later observations (Callanan et al. 1989), it was first inter- 
preted as the result of a shock caused by the interaction of 
the pulsar and companion winds. A similar situation may 
be present in PSR J2051— 0827, but the non-perfect co- 
rotation could well play an important role in deciphering 
the observed light curve. 

The variation of the projected semi-major axis of the 
pulsar can be explaind by classical spin-orbit coupling. 
Due to a significant quadrupole moment of the compan- 
ion the spin and orbital angular momenta, which are not 
aligned, couple and precess about the fixed total angu- 
lar momentum vector. The precession of the orbit results 
in the observed variation of the inclination angle whereas 
the periastron advance cannot be detected due to the ex- 
tremely small eccentricity of the system. The misalign- 
ment of the present system provides evidence that the 
neutron star received a kick at birth (cf. Kaspi et al. 1996) 
presenting further proof for the existence of asymmetric 
supernova explosions. A more detailed study will be pre- 
sented elsewhere (Lohmer et al, in prep.). 

Using timing information wc arc able to obtain a max- 
imal radius of the companion of i?cmax ^ 0.06 R^. which 
is about half the size of its Roche lobe i?L = 0.13 Rq. 
This value is consistent with that obtained by modelling 
the light curve of the companion (Stappcrs et al. 2001b). 

The quasi-cyclic variations of the orbital period show 
that there is no secular orbital decay leading to a destruc- 
tion of the binary system, which has been discussed, based 
on previous timing data, by Stappcrs et al. (1998). As 
there is no Roche lobe overflow and the timcscale to evap- 
orate the companion by the pulsar's rolativistic wind is 
~ 10^ yr (Stappcrs et al. 1998) the PSR J2051-0827 bi- 
nary system is not likely to become an isolated millisecond 
pulsar in the future. 
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